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Abstract 

The Landau-Khalatnikov time-dependent equation is applied to describe the 
crystalhzation process of the ordered vortex lattice in high temperature su- 
perconductors after a sudden application of a magnetic field. Dynamic co- 
existence of a stable ordered phase and an unstable disordered phase, with 
a sharp interface between them, is demonstrated. The transformation to the 
equilibrium ordered state proceeds by movement of this interface from the 
sample center toward its edge. The theoretical analysis dictates specific con- 
ditions for the creation of a propagating interface, and provides the time scale 
for this process. 

PACS numbers: 74.60. Ge, 64.60.My, 05.70.Fh and 64.60. Cn 
The process of formation of various equilibrium phases after a sudden change in the 
thermodynamic conditions, is a topic of wide theoretical and experimental interest. 
Obviously, the initial state created immediately after the abrupt change is a non-equilibrium, 
unstable state. The transformation of this state to the thermodynamic equilibrium state 
may proceed either homogeneously throughout the entire system, or by nucleation of a spa- 
tially localized domain of the equilibrium phase, creating a front which propagates until 
equilibrium is reached in the entire system. The latter process has been recently observed in 
the vortex system of Bi2Sr2CaCu208+5 (BSCCO) crystals. [^|| This vortex system exhibits 
a transition between an ordered and disordered phases, at a temperature-independent tran- 
sition field Bod — 400 G. [^|| High temporal resolution magneto-optical measurements [0,^ 
indicated that immediately after a sudden application of external magnetic field Ba ^ -Bod, 
a transient disordered vortex state is created, followed by a nucleation and front propaga- 
tion of an ordered vortex state. |0] The purpose of this paper is to analyze theoretically the 
crystallization process of the vortex ordered phase, i.e. the nucleation process, the creation 
of a front and its motion. 

Our analysis is based on the Landau-Khalatnikov (LK) time dependent equation |]^: 

^ = -r^^: (1) 

where \1/ and F are the order parameter and the free energy of the system, and F is the 
Landau-Khalatnikov damping coefficient. We define the order parameter of the vortex sys- 
tem in a way analogous to the definition of the order parameter in order-disorder transitions 
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in atomic solids. [§] In the latter case, the order parameter pq is a set of Fourier compo- 
nents of the atomic density taken at reciprocal lattice vectors q = G. In particular, for an 
ordered lattice phase Pg = const 7^ at g = G, whereas for a disordered state Pq = for all 
g 7^ 0. Extending this approach to the vortex order-disorder phase transition, we note that in 
small-angle neutron scattering experiments in BSCCO, ^ Bragg peaks are observed at low 
temperatures and low fields mainly in the first Brillouin zone; these peaks are smeared for 
fields larger than Bod- Thus, only one component, pc^, is sufficient to completely describe 
the order parameter, being the value of the Fourier component of the vortex density 
at the minimal vector of the reciprocal lattice. In order to describe the kinetics of the 
phase transition we allow the order parameter to be temporally and spatially dependent, 
\l/(r, t) = pG'j(r, t), assuming that \I'(r) varies slowly over the inter-vortex distance. The 
scalar real order parameter \l/(r, t) so defined, distinguishes between two thermodynamic 
solid phases of the vortex matter: \1/ = for the disordered state and \1/ = \E'o 7^ for the 
ordered state. 

In the Ginzburg-Landau formalism, the phase transition between the ordered and disor- 
dered phases may be described by a free energy density functional F: 

F = loiV^r - _ ^f3^^ + (2) 

where 7 and D are the Landau coefficients. These coefficients depend on the vortex- 
vortex and vortex-pinning interactions, and their evaluation requires a microscopic theory 
which does not yet exist. Note that Eq. 2 does not describe the whole free energy of 
the vortex system, but only that part which is varying through the phase transition, i.e. 
dependent. 

As the order- disorder vortex phase transition in BSCCO is field driven, we express the 
parameter a as a function of B: 

a = ao{l- B/B*), (3) 

where B* is a characteristic field related to the transition field Bod = B*{l+2/9p), where /i = 
ctol/P"^- Note that for a second order transition {(3 = 0), Bod = B* . For a first order phase 
transition, metastable states of the system are found between B* and B** = B*{1 + l/Ap). 
For B < B* the disordered state is unstable while the ordered state, characterized by 

^ = ^0 = ^ (1 + Vl + 4/i(l-5/5*)) , 

is stable. For B > B** the ordered state is unstable, while the disordered state with = 
is thermodynamically favorable. All the above results are deduced from the conventional 
Landau theory for phase transitions by replacing temperature with the induction B. 

In solving Eq. (|1[) we assume an initial non-equilibrium disordered vortex state = 0) 
caused by the rapid injection of the vortices through non-uniform surface barriers. |]3|JTT[| 
We show that Eq. can describe the nucleation and growth of the vortex ordered state 

= \I'o). To demonstrate this point we assume an induction distribution with a constant 
gradient |TD[ B/d, i.e. B = Ba— B{l — \x\/d), where Ba is an applied field and d is half-width 
of the sample. In this case, Eq. (|l]) can be solved analytically for both the nucleation and 
growth processes. 
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A solution for the nucleation process, i.e. the initial growth of the order parameter \1/ 
close to zero), is obtained by neglecting nonlinear terms in Eq. (Q): 



Id^f 8^^ / Ba-B B x \ ^ 

The boundary conditions dictated by symmetry is d'${x,t)/dx\x=o = 0; also we require that 
t) is a non-diverging function. The solution of Eq. (^) is then: 

oo 

^'(X, t)=Y^ Ane^-'Al{x/Xs - ^n): (5) 
n=0 

where 



A„ = Tao 



Ba-B faoY^' ( B^^"" 



B* V y \B* 



(6) 



Ai is the Airy function, c:^ = 0.685,3.9,7.06,... are the solutions of J2/?,{Sn) = J-2/ziSn) 
where Jy is the Bessel function, and Xg = lDdB*/aoB] = d[aDB*j^B\ . Here 

a£, = D'-fl f3'^d'^ is a dimensionless exchange coefficient. Note that <j„ is a constant of order 
n , growing with increasing n. 

It is evident from Eq. (|^) that only terms with A„ > play a role in the nucleation 
process. For Ba — B > B*, i.e. the entire sample is in a metastable or a stable state, all A„ are 
negative, implying that the nucleation process cannot take place. For Ba = B the induction 
at the center of the sample is zero and the rate of the nucleation process is maximum. 
Relation (|) shows that the exponent with n = yields the fastest nucleation rate, thus 
governing the nucleation process. This process may thus be approximately described by the 
first term in Eq. (^). In this approximation, the development of the order parameter during 
the nucleation process is described by the dashed lines in Fig. 1. Note that the analytical 
solution (^) describes only the first stages of the nucleation process, where the non-linear 
terms in Eq. (|l]) may be neglected. This solution ceases to apply when the value of ^ 
approaches \l/o, i.e. after a time period of order I/Aq. The width of the ordered domain is 
then If ~ Xs (1 -|- <^o) ~ a;^. The condition for appearance of localized domain in the sample 
center may be then obtained from the inequality Xg << d or: 

i- » -• 

B* fx ^ ^ 

If this condition is not satisfied, then homogeneous transformation of the unstable phase takes 
place. Otherwise, a sharp front will develop separating between the nucleating ordered phase 
and the initial unstable disordered phase, as described above. Thus, when the induction 
gradient is large enough, we expect the appearance of a sharp interface between the growing 
stable (ordered) phase and the retreating unstable (disordered) phase. 

In describing the growth process, i.e. the movement of the interface between the ordered 
and disordered phases, non linear terms in Eq. |l| must be taken into account. We express the 
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linearly varying function B{x) as: B = B f + {B / d){x — x f) , where Bf = — B) + XfB/d is 
the induction at the front located at Xf. In this notation: a = 1 — + {B / d){x ~ x f)]/ B* . 

Eq. (|l]) can be written in the reference frame of an observer moving with the front, by 

t 

introducing a new variable ^ = x — Xf (t) and defining Xf (t) = xq + f Vf{t')dt', where Vf 



is the time-dependent front velocity and Xq is a constant. With the new set of independent 
variables Bf) Eq. (|l]) becomes: 



9^ Bd^^X ^92* / Bf{xf{t)} B^ . . 
One can solve this equation analytically provided the front width A << dBf/ B. In this case, 



the terms {vf/T){B/d){d^/dBf) and -ao{B^/ dB*)^ , may be neglected, The solution 
of Eq. d) is then [T|]: 

vl> = vl>o(l + expfiH , (9) 



A 



where 



A^ 

A2 = ° (10) 

2/i(l - Bf/B*) + 1 + + 4/i(l - 

is the front width. The front velocity Vj = dxf/dt is: 

^/ = ^0 / (11) 



Here 



^^0 



2/i(l - + 1 + - Bf/B*) 



T^/WWf = Al = AD^/f3'' = Ad'^ao, (12) 



and Bf = Bf {xf{t)}. The solid lines in Fig. 1 show \1/ given in Eq. (|^) for different locations 
of the front Xf, describing the propagation of the ordered phase. 

We turn now to discuss the front velocity Vf, Eq. (|TI|), and the front width A, Eq. 
We first note that Vf and A do not depend explicitly on time or applied field but on Bf, 
the local induction at the front. Several important conclusions may be drawn from these 
equations: 

1) As Bf approaches Bod the velocity is approaching zero (the Vf{Bf) dependence is 
described by the sohd line in Fig. 2). 

2) The motion of the front toward the sample edge is accompanied by an increase of the 
induction Bf at the front, resulting in a decrease of the velocity with time. 

3) The front width A decreases with the increase of (3, implying that for a 'stronger' 
first order transition the front is steeper. Also from Eq. (|TD[) it is obvious that the exchange 
coefficient D causes the front to be smeared. In addition, increasing D and/or the damping 
coefficient F, results in an acceleration of the front motion (see Eq. |Tl|). 
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So far we have demonstrated dynamic coexistence of ordered and transient disordered 
vortex phases, with a sharp interface between them, assuming time-independent induction 
distribution with a constant gradient. In high-temperature superconductors, however, the 
induction distribution varies significantly with time due to flux creep. In addition, one may 
expect different flux creep laws for the different vortex phases. As a result, the nucleation 
and growth of the ordered vortex phase are manifested experimentally by the appearance 
of a break in the induction profile and movement of this break toward the sample edge. |0] 
The location of the break is expected to coincide with the location of the moving front of 
the order parameter. 

To demonstrate this scenario we solved numerically the LK equation (|I|), allowing for 
flux creep. We assume a \l/-dependent local current density 

ji^,t) = Mt)(i-§-^+Mt)§-, (13) 

where Ji{t) and J2{t) are the current densities in the disordered and in the ordered phases, 
respectively. |]14| Based on the experimental observations we further assume Ji{t) = 



Joi(l + t/ti) and J2{t) = ^02(1 + t/t2) The induction B (and thus the coefficient a) 
may now be expressed in terms of the order parameter by using Maxwell equation: 

d 

B{x,t) = Ba — Atx/c J J{'^{y,t),t)dy. As before, we assume an initial disordered phase 

X 

throughout the entire sample, = 0,x) = 0. 

In order to solve equation ([^) numerically we define dimensionless parameters: b = B/B*, 
J = 47TJd/{cB*),x' = x/d, t' = t(3^T/-f, ^' = ^/^o{B**) = 2^7//?. Eq. (|ID then becomes: 

^ = + ^(1 - bW))^' + - + /'(x', t'), (14) 

where /' = 2/7/(/?ao) is a dimensionless noise which must be introduced in the numerical 
solution. 

The values of the (dimensionless) parameters used in the numerical calculations are based 
on experimental measurements. In particular, from the fit of Eq. (|11]) to the experimental 
data of Vf{Bf) in Fig. 2 we estimate = 1.5, thus Bod/B* = 1.148, B**/B* = 1.166. A 
value of 10""^ for ao = (2/ifo/(rao(i))^ (see Eq. is estimated from the experimental value 
Vq ~ 20 ixm/s obtained from the same fit; a value of Aq ~ Ta^ ~ 10 s^^, is estimated from 
the time elapse between switching-on of the external field and the appearance of a break in 
the induction profile. A typical value of ci ~ 300 iim for the sample half width was used. 
Based on the analysis of magnetic relaxation we take ai = 0.3, 02 = 0.5. In addition, a 
noise level of f^^^ = 10~^ is assumed. The system of equations completed by boundary and 
initial conditions has been solved numerically utilizing the Euler method. The unit space 
interval was divided into 200 segments, and time step of 2.5 ■ 10~^ (in dimensionless units) 
was used providing a stability of the numerical procedure. The results for Ba/ B* = 1.1 are 
shown in Fig. 3. The upper case shows the spatial dependence of the order parameter at 
different times. The nucleation appears at the sample center at t ~ 10, forming a sharp 
front which propagates toward the sample edge. |T^ The lower case of Fig. 3 shows the time 



evolution of the induction profiles during the nucleation and growth processes. A sharp 
break in the profiles appears at the location of the front of the order parameter after the 
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nucleation is completed. As expected, the break in the induction profile and the front of the 
order parameter move together toward the sample edge. Note that a break in the induction 
profile can be observed outside the region of phase metastability (i.e., for Bf < B*). 

The theoretical predictions described above are confirmed experimentally in BSCCO 
crystals. In particular, breaks in the induction profiles were recorded following a sudden 
application of external field of intensity close to Bod- This break moves toward the sample 
edge at a velocity which depends only on Bf, the value of the induction at the break. Thus, 
the dependence of on 5^ is not affected by magnetic relaxation. As shown in Fig. 2, the 
analytical curve (solid line) of Vf{Bf), Eq. (|TTp, is in good agreement with the experimental 
results (open symbols). Moreover, numerical results (solid symbols in Fig. 2) for Vf{Bf) for 
different applied fields also show a good agreement with the analytical curve demonstrating 
that magnetic relaxation does not affect the dependence of Vf on Bf. 

Finally, we note that two velocities govern the vortex dynamics in the process of the 
phase transformation: The interface velocity Vf, and the flux velocity, vp, due to creep. The 
effect of the latter on the shape of the interface must be taken into account in close vicinity 
of Bod where Vf — > 0. 
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Figure captions 



Figure 1. Nucleation and growth of the order parameter. The nucleation process is 
demonstrated by the dashed curves calculated from Eq. (|) for An = ^oi^n.o at times 
Aot = 8.12, 9.86, 11.02. The solid lines, describing the growth process, are calculated from 
Eq. at different locations Xf/d = 0.3, 0.5, 0.7, 0.9. 

Figure 2. Experimental data (open symbols) of Vf{Bf) for different applied fields, taken 
from Ref. 0, a fit (solid line) to Eq. |Tl| and results of numerical solutions (solid symbols) 
for the indicated applied fields . The fit was done with two fitting parameters vq and /i, after 
Bod was estimated to be approximately 400 G, corresponding to an induction value where 
the velocity is going to zero. 

Figure 3. Order parameter (upper case) and induction profiles (lower case) for ji = 
1.7(1 + t')~^'^ and j2 = 1-5(1 +t')~^'^. The profiles are shown for dimensionless times t' = 1, 
4, 9, 10, 11, 12, 14, 16, 18, 22, 27, 35, 44, 54, 65, 79, 99. 
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